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FOREWORD 


Recovery  of  payloads  from  high  altitude  often  requires  the 
deployment  of  parachutes  on  trajectories  that  are  essentially 
vertical  and  toward  the  earth.  The  parachute  opening  shock  force 
developed  by  a  parachute  deployed  in  this  manner  exceeds  the 
opering  shock  force  of  the  same  parachute  system  deployed 
horizontally  at  the  same  altitude  and  velocity.  As  the  deployment 
trajectory  angle  varies  from  horizontal  to  vertical,  the  opening 
shock  force  increases  to  a  maximum.  Hence,  the  vertical 
trajectory  has  special  significance.  This  report  develops  an 
analysis  to  permit  the  calculation  of  opening  shock  forces  in 
vertical  deployment.  The  validity  of  the  rule  of  thumb  that  "the 
vertical  deployment  opening  shock  force  exceeds  the  horizontal 
shock  force  by  one  "g"  is  tested,  and  criteria  are  developed. 

Examples  are  used  to  demonstrate  applications  of  the  approach 
and  a  method  of  calculating  the  inflation  time  of  solid  cloth 
parachutes  in  vertical  fall  is  developed. 

The  author  wishes  to  express  his  appreciation  to  Mr.  Hensel 
Brown  of  the  Strategic  Systems  Department  for  his  valuable 
assistance.  Mr.  Brown  developed  the  computer  programs  and 
conducted  numerous  background  calculations  for  this  report. 
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Dr.'-O.  E.  GOELLER,  Head 
Underwater  Weapons  Division 
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INTRODUCTION 


Parachutes  are  used  to  reduce  the  impact  velocity  of  many 
falling  objects.  Often  the  parachute  systems  are  initiated  while 
moving  in  a  trajectory  which  is  parallel  to  the  earth's  surface. 
Many  methods  of  calculating  parachute  opening  shock  force  use  the 
horizontal  attitude  as  a  basis  for  analysis.  However,  in  the 
recovery  of  systems  which  have  been  allowed  to  free  fall  from  high 
altitudes,  the  trajectory  deviates  markedly  from  the  horizontal 
and  approaches  the  vertical  direction.  The  opening  shock  force  of 
a  vertically  towa rd-the-ear th  deployed  system  is  known  to  exceed 
the  horizontal  deployment  shock  force  for  the  identical  system  and 
operational  altitude  and  velocity. 

Vertical  trajectories  are  of  special  interest  since  the 
opening  shock  force  is  a  maximum  value.  The  purpose  of  this 
report  is  to  develop  a  method  for  calculating  the  drag  area, 
velocity  and  shock  factor  profiles,  opening  shock  force,  time  of 
occurrence  of  the  shock  force  during  canopy  inflation,  and  the 
distance  of  fall  for  parachutes  deployed  in  a  vertical  toward-the- 
earth  trajectory.  This  report  is  an  extension  of  Reference  1, 

NSWC  TR86-142,  "Notes  on  a  Generic  Parachute  Opening  Force 
Analysis",  which  describes  a  technique  of  calculating  the  opening 
shock  force  of  several  types  of  horizontally  deployed  parachutes. 
In  the  interest  of  clarity  readers  should  familiarize  themselves 
with  Reference  1,  as  it  is  to  be  used  extensively  in  this  report. 

A  method  of  calculating  the  inflation  time  of  solid  cloth 
parachutes  in  vertical  fall  is  developed  and  horizontal  and 
vertical  inflation  times  for  identical  deployment  conditions  are 
compared.  Criteria  required  for  inflation  time  analysis  of  other 
types  of  parachutes  are  presented. 

There  is  a  rule  of  thumb  used  in  calculating  the  opening 
shock  force  of  vertically  deployed  parachutes.  The  rule  states 
that  "the  vertical  mode  opening  shock  force  is  equal  to  the 
horizontal  mode  opening  shock  force,  plus  one  g".  Examples 
demonstrate  that  the  vertical  deployment  mode  inflation  reference 
time  is  less  than  the  horizontal  deployment  mode  inflation 
reference  time  for  the  same  system,  velocity  and  altitude.  The 
theory  indicates  that  the  vertical  deployment  mode  opening  shock 
forces  exceed  the  rule  of  thumb. 
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APPROACH  i 
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The  analysis  is  based  upon  the  application  of  Newton's  second 
law  of  motion  to  particle  trajectories  of  inflating  parachutes  in 
the  horizontal  and  toward-the-earth  vertical  modes.  Equations  for 
the  horizontal  deployment  of  parachutes  with  dynamic  drag-area 
signatures  of  the  form  CDS/CDS0  = (1-r)  (t/t0)D  +7  were  developed 
in  Reference  1.  The  same  dynamic  drag-area  signatures  are  to  be 
incorporated  into  the  vertical  deployment  analysis.  The  exponent 
j  determines  the  type  of  parachute.  A  value  of  j=l  is  indicative 
of  geometrically  porous  ribbon  and  ring  slot  types  of  canopies. 
Solid  cloth  parachutes  have  a  value  of  j=6. 
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FIGURE  1.  VERTICAL  TRAJECTORY  POINT  MASS  FORCE  SYSTEM 
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DEVELOPMENT  OF  THE  VELOCITY  -  TIME  EQUATIONS 
With  reference  to  Figure  1. 


dv 

CoS 

CdSo 

dV 


=,,-r,(i:)'+r 


(2) 


where  the 
1  to  be  a 


quantity  2W/(pg  Vs  t  Cq  S  )  has  been  shown  in  Reference 
ballistic  mass  ratio  scale°factor. 


Let  M  = 


2W 

pg  V*t0CoSo 


(3) 


h =1  A* -*&:(" -T]  iz)‘ "h* 


V  =  V,  t  =  0 


(4) 


DEVELOPMENT  OF  THE  DISTANCE-TIME  EQUATION. 

Application  of  the  variable  dynamic  drag  area  and  ballistic  mass 
ratio  equations  to  the  velocity  distance  form  of  equation  (1) 
yields. 


s  v 

Ids  =  I 

S  —  0  v  =  v, 


L _ 

\  MV, t0 


VdV _ 

("-T'id'  +  T)) 


(5) 


At  any  instant  during  the  inflation  process  the  shock  factor 
lias  been  shown  to  be: 


X  ;  - 


CoS 

CdS0 


(6) 
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and  the  instantaneous  shock  force  is: 


Fi  =XiF, 


(7) 


where  Fa  is  the  steady-state  drag  force  of  the  fully  open 
parachute  at  the  line  stretch  velocity,  Vs . 


F»  =  CoSo  (8) 


The  opening  shock  performance  in  the  vertical  and  horizontal 
deployment  modes  shall  be  demonstrated  by  examples. 

Example  1:  Compare  the  horizontal  and  vertical  deployment 
mode  opening  shock  forces  for  a  flat  circular  solid  cloth 
parachute  of  35  feet  DQ  diameter  and  30  gores.  Deployment 
conditions  are: 


a . 

System  weight  =  200  lb. 

b. 

Density  altitude  =  3,000 

f  t . 

c . 

Line  stretch  velocity  = 

340 

f  ps . 

d . 

Drag  coefficient  =  0.75 

e . 

Dynamic  drag  area  exponent, 

3  =  6 

f . 

Initial  drag  area,  r=0 

HORIZONTAL  OPENING  SHOCK  FORCE 

a.  Parachute  surface  area 

_  ir  _  2 
So  —  q  Do 

So  =  y <35l2 

So  -  962.1  ft2 
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b.  Parachute  drag  area: 

CdSo  =  0  75  x  962.1 
CdSo  =  721.6  ft2 


c . 

d. 


Air  density,  p  =  0.0021753  slugs/ft' 
Steady-state  drag  force 

Fs  =  y/>VsCDSo 

*j 

Fs  =  —  X  0.0021753  (340)2  721.6 
F$  =  90726.71b 


e. 


In  order  to  determine  the  ballistic  mass  ratio,  the 
inflation  time  tQ  must  be  known.  Two  methods  for  calculation  of 
the  horizontal  deployment  reference  time  are  available.  One  is 
via  equation  (9)  from  page  A-7  of  Appendix  A  and  the  other  via 
equation  (10)  from  page  5-61  of  Reference  2. 


to  — 


2W 


pgVjCoSo 


gpVo 

CdSo 

2W 

am„  -a„  k/C  P  ',V! 

- 

\  2  / 

1 

< 

_ 

- 

(9) 


From  Table  2,  page  A-15  of  Appendix  A,  the  following  inflated 
shape  data  were  obtained  for  a  30-gore  solid  cloth  flat  circular 
parachute.  See  Figure  24  of  Appendix  A  for  shape  nomenclature. 
Steady-state  inflated  canopy  radius: 


~~  -  0.668:  ™  =  0  827.  £  =  0  6214;  t 
D0  a  a  a 

0.668  x  35 


5  =  11  69  ft 


0  7806 
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fi 

F 


C'. 


Steady-state  mouth  area: 

_,r  /N/a  —  b/SN2a 

AMo  =  ™*[l-(-  —  )  J 

/  0.827  -  0.6214 \2 


AMo  =  -n-(11.69l2[l 
Amo  =  399.53  ft2 


0.7806 


] 


Steady-state  canopy  volume  of  air  to  be  collected,  geometric. 

,,  2  - 3 r b  ,  b'“] 

Vo  =  Jtra  [T+1J 

Yo=  j  m  (11. 69) 3  [o. 6214  +  0.7807] 


Yo  =  4690.83  ft3 


The  MIL-C-7020,  type  III  canopy  cloth  airflow  constant  is  k=  1.46 
and  the  average  pressure  coefficient  is  taken  to  be  C.P.  =  1.7 
then: 


to  — 


14x200 


K1  = 


0.0021753  x  32.2  x  340  x  721  58 
32.2  x  0.0021753  x  4690  33 


[•"  -’] 


2x200 


721.58 


399  53  -  962.1  (1  46) 


^  0.0021753  x  1.7  y  2 


t0  =  0  773  sac . 


Knacke,  Reference  2,  presents  the  inflation  time  as 

r.Do 


tf  - 


t» 


V, 

8  x  35 
340 


tt  =  0.824  sec. 


where  n=8  for  solid  cloth 
parachutes 
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While  equation  (10)  is  more  convenient  than  (9)  to  use,  equation 
(9)  gives  the  effects  of  altitude,  mass,  and  cloth  rate  of 
airflow;  both  values  are  estimates  of  the  inflation  time. 

Equation  (9)  was  theoretically  derived.  Equation  (10)  was 
empirically  developed  by  other  experimenters  from  numerous  field 
test  data.  The  parachute  system  of  Example  1  for  standard  air 
density  at  an  altitude  of  1000  feet  yields  values  of  to=0.839 
seconds  and  tf=0.824  seconds.  The  closeness  of  the  theoretical  tQ 
and  the  empirical  tf  values  indicates  that  the  application  of  the 
criterion  listed  below  is  a  reasonable  approach  to  the  inflation 
time  calculation. 

a.  Average  pressure  coefficient 

b.  Cloth  permeability 

c.  The  trajectory  velocity  is  taken  as  the  canopy  inflow 
velocity 


d.  The  ratio  of  the  instantaneous  mouth  area  to  the  steady- 
state  mouth  area  is  equal  to  the  instantaneous  drag-area  ratio 

e.  The  ratio  of  the  instantaneous  pressurized  canopy  area  to 
the  canopy  surface  area,  S0,  is  equal  to  the  instantaneous  drag- 
area  ratio 

f.  The  volume  of  air  to  be  collected  during  the  inflation 
process  is  represented  by  .V  . 

As  a  demonstration  of  the  cloth  rate  of  airflow  effects  in 
equation  (9),  replacement  of  the  canopy  cloth  of  example  1  with  a 
high  rate  of  airflow,  three  momme  silk  cloth  (K=7.43)  yields  an 
inflation  reference  time  of  t0=99.76  sec.  under  the  example  1 
operational  conditions.  This  time  exceeds  the  time  of  fall  from 
3,000  feet  and  the  system  would  appear  to  an  observer  as  a  partial 
inflation.  The  same  system  deployed  horizontally  at  20,000  feet, 
at  constant  dynamic  pressure,  has  an  inflation  reference  time  of 
t0=1.5  sec.  When  the  area  term 


Amo  -  Asok 


=  O 


equals  zero,  the  inflation  reference  time  is  infinite,  (eq.  9). 

As  the  area  term  approaches  zero,  t0  increases  exponentially.  It 
is  possible  for  the  inflation  time  to  be  extended  to  the  point 
that  the  parachute  system  is  not  practical  even  though  the  area 
ratio  is  not  zero.  Parachute  canopies  without  cloth  airflow,  K.=  0, 
always  inflate. 
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When  using  generalized  expressions  such  as  equation  (10),  it 
is  important  to  know  the  mode  of  operation  for  which  the  formula 
was  determined.  The  identical  parachute  used  in  finite  mass, 
intermediate  mass,  or  infinite  mass  systems  experiences  a 
variation  in  inflation  time  due  to  its  mode  of  use.  At  the  time 
tQ  all  inflations  have  collected  the  same  volume  of  air.  Finite 
mass  inflations  experience  a  significant  velocity  reduction  during 
the  inflation  process  which  limits  the  rate  of  airflow  into  the 
canopy  and  extends  the  inflation  time.  Infinite  mass  deployments, 
on  the  other  hand,  can  be  nearly  constant  velocity  which  maintains 
a  higher  rate  of  mass  flow  into  the  canopy  and  decreases  the 
inflation  time.  Actual  constant  velocity  infinite  mass  operation 
produces  the  minimum  inflation  time.  The  inflation  reference 
times,  for  the  parachute  of  Example  1,  are  plotted  in  Figure  2  as 
a  function  of  weight-to-drag-area  ratio  for  equation  (9).  The 
finite  mass  and  infinite  mass  limiting  ballistic  mass  ratios  for 
solid  cloth  parachutes  are  referenced  on  the  graph. 

The  trends  of  Figure  2  show  that  the  principal  variation  of 
the  inflation  reference  time  occurs  in  the  finite  mass  mode  of 
operation.  By  the  transition  from  finite  mass  to  low  intermediate 
mass  a  near  minimum  inflation  reference  time  is  reached .  This 
quickly  converges  to  an  essentially  constant  value  for  the  high 
end  of  the  intermediate  mass  and  infinite  mass  regimes.  Similar 
trends  for  other  types  of  parachutes  are  a  reasonable  extension  of 
the  premise  since  it  is  based  on  variable  rates  ot  flow  filling  a 
constant  volume. 

Ballistic  Mass  Ratio  for  tQ  =  0.773  sec. 


M  = 

M  - 


_ 2W 

pgVjtoCoSo 

2  x  200 _ 

0.0021753  x  32.2  x  340  x  0.773  x  721.58 


M  =  0.0301 


Since  the  ballistic  mass  ratio  is  less  than  the  limiting 
value  of  0.1907,  the  inflation  performance  is  in  the  finite  mass 
range  of  operation  and  tQ  =  tf.  The  following  formulae  were 
obtained  from  Table  12  ot  Reference  1  for  values  of  j=6  and 
r  *  0  . 


1  0 
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The  time  of  occurrence  of  the  maximum  shock  force  during 
canopy  inflation. 


t  ^  /  21M  y 

to  01  X'rn.*  "  \  4  / 

t  _  /  2 1  X  0.0301  V  ’ 

xl  =  - 3 - ) 


—  at  X,  =  0.7685 

t0  n‘<  m»« 


The  maximum  shock  factor: 

y  =  21  f21M\67 
49  \  4  / 

v  16  /  21  x  0  0301  \6  7 

Xim.,  ~  4§  \  4  ) 

Xlnl„-  0.0672 

Maximum  shock  force  in  the  horizontal  mode  of  deployment. 

F  c  Y  F 
X'm«  » 

F„,„  *  0.0672  x  90,726.7 
mix 

" 61  ' 

max 


The  magnitude  of  th 
inflation  reference  time 
regime  of  operation.  In 
the  measured  inflation  t 
which  results  in  a  varia 
horizontal  opening  shock 
technique  for  assigned  v 
results  are  plotted  in  F 
forces  and  the  "lg  added 
tQ  and  tf  for  example  (1 
performance. 


e  opening  shock 
of  the  parachut 
the  field  there 
ime  from  the  nom 
tion  in  opening 
forces  was  made 
a  ^s  of  t0  from 
igure  3  for  the 
"  rule  of  thumb. 
)  mark  the  range 


force  varies  with  the 
e  in  the  finite  mass 
i'  often  a  variation  in 
inal  calculated  value, 
shock  force.  A  survey  of 
via  the  foregoing 
0.5  through  1.0.  The 
horizontal  deployment  mode 
The  particular  values  of 
of  the  nomina]  expected 


t  1. 1  J1.1  ii- 
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A  method  of  calculating  "t0"  in  the  vertical  mode  of 
deployment  is  required  so  that  the  horizontal  and  vertical 
inflation  reference  times  can  be  compared  for  identical  systems 
deployed  at  the  same  altitude  and  velocity. 

The  inflation  reference  time  '*t0“  in  the  vertical  deployment 
mode  can  be  determined  from  the  basic  mass  flow  equation  used  to 
calculate  the  horizontal  deployment  "tp" .  Following  the  method  of 
Section  IV  of  Appendix  A  for  solid  cloth  parachutes. 


*  rnmflow  ~  mourflow 


VA  -  fi  AtP 
dt 


Instantaneous  canopy  mouth  area. 


Am  -  Amo 


FOR  : =  O 


Instantaneous  pressurized  canopy  cloth  area. 


FOR  r  =  o 


Instantaneous  cloth  rate  of  airflow. 


p*k/C-(>^  v2n  c.f.s./ft2 

\  2  ) 


the  instantaneous  velocity  is  determined  by  equation  (4)  for 


v 

V  =  Idv 
V  =  Vs 


(4) 


Where 
T=  0. 
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The  calculated  canopy  volume,  calc'  is  determined  from  equation 
(11)  . 


A  program  for  calculating  tQ  for  solid  cloth  parachutes  and 
the  opening  shock  force  profile  during  the  inflation  of  several 
parachute  types  is  provided  in  Table  1.  Equations  (4),  (11),  and 
(12)  are  programmed  together  with  the  vertical  deployment  opening 
shock  equations  (2)  through  (8)  in  FORTRAN  IV  language.  The 
included  examples  were  calculated  via  the  program  using  a  VAX  780 
computer.  The  program  operates  in  two  modes.  Mode  1,  for  solid 
cloth  parachutes,  calculates  the  vertical  deployment  reference 
time  tc  for  the  parachute  system  parameters  and  operational 
deployment  data,  and  then  calculates  the  opening  shock  profile 
during  inflation.  A  typical  data  print  out  is  shown  in  Table  2. 

It  is  necessary  to  estimate  an  initial  value  of  ”t0"«  The  program 
calculates  the  canopy  volume  for  the  estimated  time  and  compares 
the  ¥0  calc  to  the  volume  derived  from  the  canopy  geometry.  If  the 
calculated  volume  is  not  within  specified  limits,  the  program 
adjusts  ‘’to”  by  equation  (12)  and  reiterates  the  program  until  the 
calculated  volume  is  within  the  specified  limits. 


Mode  2  of  the  program  calculates  opening  shock  profiles  for  input 
values  of  t0.  Mode  2  analysis  of  other  types  of  parachutes  is 
possible  by  the  selection  of  the  proper  values  of  ,‘jM  (1/2,  1,  2, 

3,  4,  5,  or  6)  and  "r".  The  opening  shock  force  variation  for 
examples  (1)  and  (2)  are  plotted  in  FigureE  3  and  4.  The  nominal 
t0  for  n=0 . 632  was  calculated  by  the  program  in  mode  1  and  the 
force-time  survey  was  calculated  in  mode  2. 

Figure  3  illustrates  that  parachutes  deployed  in  a  vertical 
toward-the-earth  trajectory  inflate  faster  than  the  same  system 
deployed  horizontally  at  the  same  altitude  and  velocity. 

Inflation  reference  times  for  parachute  types  other  than 
solid  cloth  canopies  can  be  developed  from  the  mass  flow  equation. 
This  requires  that  the  flow  through  the  canopy  be  expressed  in  a 
form  similar  to  the  solid  cloth  canopy  cloth  permeability,  P, 
where  the  rate  of  flow  per  unit  area  is  a  function  of  the  pressure 
differential  across  the  cloth  or  grid. 

!  i 
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TABLE  1.  INSTANTANEOUS  ORAG  AREA,  VELOCITY,  OPENING  SHOCK  FORCE.  AND 
OISTANCE  OF  FALL  OF  A  PARACHUTE  DEPLOYED  IN  VERTICAL  FALL. 


THIS  PROGRAM  CALCULATES  THE  INSTANTANEOUS  DRAG  AREA.  VELOCITY, 

OPENING  SHOCK  FORCE,  AND  DISTANCE  OF  FALL  OF  A  PARACHUTE  DEPLOYED  IN 
VERTICAL  FALL 

THE  PROGRAM  OPERATES  IN  TWO  MODES 

MODE  1  -  CALCULATES  THE  INFLATION  TIME  AND  PERFORMANCE  PROFILES  FOR  SOLID  CLOTH 

PARACHUTES  (TO  INPUT  AS  INITIAL  ESTIMATE)  (IOPT  =  1) 

MODE  2  -  CALCULATES  THE  PERFORMANCE  PROFILES  FOR  VARIOUS  TYPES  OF  PARACHUTES  (J) 

INFLATION  TIME  INPUT  IS  REQUIRED  (IOPT  =  2) 


INPUT  IOPT 


1  (FOR  MODE  1 ) 

2  (FOR  MODE  2) 


INPUT  NEEDED  FOR  BOTH  MODES 


RHO 

VS 

CDSO 

TO 

U 

J 


AIR  DENSITY  AT  GIVEN  ALTITUDE  (SLUGS/FT3) 

VELOCITY  AT  SUSPENSION  LINE  STRETCH  (FT/SEC) 

DESIGN  DRAG  AREA  (FT2) 

1  OPT  *  1  INITIAL  CUESS  FOR  INFLATION  REF  TIME  (SEC) 
I  OPT =2  ACTUAL  INFLATION  REFERENCE  TIME  (SEC) 
UEIGHT  (LHS) 

=6  FOR  FLAT  CIRCULAR  PARACHUTE 
=1  FOR  RIBBON  TYPE  OF  PARACHUTE 


INPUT  NEEDED  FOR  IOPT  =  1  ONLY 


AMO  -  STEADY-STATE  MOUTH  AREA  ( FT2 ) 

ASO  -  CANOPY  DESIGN  SURFACE  AREA  (FT2) 
K  -  CLOTH  PERMEABILITY  CONSTANT 

CP  -  PRESSURE  COEFFICIENT 

N  -  CLOTH  PERMEABILITY  EXPONENT 

VO  -  CEOMETRIC  VOLUME  (FT3> 


REAL«4  N 
T  ODEN*- 1 OOOOO 
5  PRINT  *,  ' INPUT  IOPT  ' 

READ(5.  *,  END*-100)  IOPT 
PRINT  *,  'INPUT  RHO,  VS,  CDSO,  TO,  W,  J  ' 
READ ( 5,  *  >RHO,  VS,  CDSO,  TO,  W,  J 
IF ( IOPT  EQ  21  GO  TO  3 
PRINT  *,  'INPUT  AMO,  ASO.  XK,  CP,  N. VO  ' 
READ (  5,  *  )  AMO,  ASO,  XK,  CP,  VO 
3  DT=TO/TODCN 
TAU»0 
C-  32  2 
S=0 
X*=TAU 

CDS=TAU*CDSO 

FS=  5*RHQ*VS**2 *C DSO 

f«tau*fs 

V0L*0 
I P ASS=0 

1  u 
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TABLE  1.  (CONT.) 


DT=  10  *DT 

IF ( I OPT  EQ  2)  GO  TO  8 

7  V=VS 
V0L*0 
1-0 

XM=<2  *W)/<RHO«C*CDSO*VS*TO> 

6  DV=(G-V**2/<  XM*VS*T0)*<T/T0)**J)*DT 

DVOL=<  V#AMO* ( T/TO ) ** J-ASO* ( T / TO  > **J*XK* ( CP*RH0/2  >**N 
1  *V**<2  *N) ) *DT 
VOL=VOL+DVOL 
V-V+DV 
T*T+DT 

IF  (  T  CE  TO)  GO  TO  4 
GO  TO  6 

4  IF (VOL.  GT.  (VO-10  )  AND  VOL  LT  (V0+10  >)  GO  TO  6 
T0*T0*< VO/VOL  > 

IF< IPASS  GT  50)  STOP 
IPASS=IPASS+1 
GO  TO  7 

8  CONTINUE 
T =0 

v*vs 

DT =DT / 10 

WRITE<6,  10)  RHO,  VS.  CDSO,  W,  TO.  XM,  J,  FS 
10  FORMATOHl,  'DENSITY* ' 1PE12  5,  2X.  'V(S)*'0PF8  1.2X,  'CDSO* 'F8  2, 

1  2X,  'W*'F7  1 . 2X,  'T(0)='F6 .3,  2X.  'M='1PE12  5.2X,  'U='I2,2X, 

2  'FS* ' 1PE12  5) 

IF ( I OPT  EQ  1 )  WRITE<6.  11)  AMO,  ASO,  XK, CP, N, VO, VOL 
1  1  FORMAT  <  1H0,  4X,  'AMO  *  '  ,  0PF8  2.  7X,  'ASO*  0PF8  2,  5X,  'K*  '  ,  0PF8  3.  IX, 

4  'CP*  ' , 0PF7  3.  3X,  'N=  OPFA  3.  IX,  'VO* ',  1PE12  5,  7X,  'VOL* ' 1PE12  5) 

WR I TE  <  6, 20) 

20  FORMAT < 1H0, 2X ,  'TIME '2.  ,  'TIME  RATIO',  2X,  'VEL  RATIO', 2X,  'DRAG  AREA 
1RAT 10 ' , 2X , 'SHOCK  FACTOR', 3X,  ' SHOCK  FORCE'. 2X, 'DIST  OF  FALL') 

30  FORMAT (IX, F6  2, 2F 1 2  5.F17  7 ,  1 P2E 1 4  5,  OPF 1 5  2 ) 

LCOUNT-5 

33  IF ( T/TO  CT  1  51)  GO  TO  5 
I F  <  LCOUNT  LT  54)  GO  TO  38 
WR ITE<  6, 37) 

37  FORMAT ( 1H1 ) 

WRITE (6. 20) 

LC0UNT=3 

38  LCOUNT =LCOUNT  + 1 

WR  I  TE  (  6,  30)  T,  T/TO,  V/VS,  CDS/CDSO,  X,  F,  S 
TOLD*T ♦  02 

IFIDT  6T  0  02)  TOL D*=T +DT 
40  T*T  +DT 

IFU  GT  TOLD)  T  =  TOLD 

TCVTO=  < 1  -TAU)*(T/TO)**J  +  TAU 

CDS*CDS0*T0VT  0 

DV*(C-V**2*T0VT0/ <  XM*VS*TO)  )*DT 
DS*V*DV/ <C-V**2#T0VT0/ ( XM*VS*TO ) ) 

V*V*DV 

S*S+DS 

X*< V/VS >**2»CDS /CDSO 
F*  X*FS 

IF  ( T  LT  TOLD)  GO  TO  40 
GO  TO  35 
100  STOP 
END 

15 
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The  various  t0  and  tf  times  noted  on  the  graph  are  the 
particular  solutions  for  the  deployment  conditions  of  example  1. 
All  of  the  ballistic  mass  ratios,  M,  in  the  calculations  are  less 
than  the  limiting  value  (Ml=0.1907)  for  finite  mass  operation. 

Also,  the  calculated  opening  shock  forces  are  less  than  the 
suspension  line  constructed  strength  (Fc=16,500  lb)  which 
indicates  a  practical  system  for  the  opening  shock  requirement. 
Examination  of  the  inflation  times  for  equations  (9)  and  (11)  for 
n=0.5  demonstrates  that  the  vertical  deployment  parachute 
inflation  time  is  reduced  even  though  the  system  parameters  of 
deployment  altitude  and  velocity  were  identical.  The  reason  for 
this  is  that  parachutes  are  energy  dissipators.  In  the  horizontal 
deployment  mode  the  additional  added  potential  energy  of  the 
system  weight  falling  a  short  distance  during  inflation  is  not 
included  in  the  calculations.  In  the  vertical  mode  considerably 
more  potential  energy  is  absorbed  into  the  system  resulting  in 
less  trajectory  velocity  reduction  during  inflation  and  hence  a 
shorter  inflation  reference  time. 

In  field  tests  the  inflation  reference  time  often  varies  from 
the  calculated  value.  One  cause  of  this  is  the  variation  of 
airflow  through  the  cloth  as  it  is  -manufactured  and/or  handled. 

The  cloth  airflow  constants  (k,  n)  would  have  to  be  true  average 
values  for  each  parachute  which  requires  so  many  fabric  tests  as 
to  be  unreasonable.  The  effects  of  varying  n  for  a  constant  k  are 
shown  in  Figure  3  for  equation  (11).  As  n  increases  from  0.5  the 
rate  of  cloth  airflow  is  augmented,  extending  the  inflation 
reference  time  and  reducing  the  opening  shock  force. 

For  any  constant  inflation  reference  time,  the  vertical 
deployment  mode  opening  shock  force  exceeds  the  "rule  of  thumb" 
value. 

Example  2:  Investigate  the  feasibility  of  retarding  a  2000- 
lb.  payload  using  the  parachute  system  of  example  1. 

Figure  4  presents  a  horizontal  deployment  analysis  using  the 
methods  of  example  1  for  assigned  values  of  tQ  from  0.3  to  0.8 
seconds.  The  vertical  deployment  analysis,  using  the  program 
modes  of  Table  1  was  conducted  using  the  same  tc  range.  The 
particular  inflation  reference  times  for  the  conditions  of  example 
2  was  computed  by  equation  (9)  and  equation  (11).  The  increased 
system  weight  caused  a  reduction  in  all  values  of  t0  as  expected 
from  Figure  2.  The  inflation  time  calculated  from  equation  (10) 
is  unchanged  from  Figure  3  since  equation  10  does  not  adjust  for 
weight  variation. 
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The  analysis  balJistic  mass  ratios  exceed  the  finite  mass 
limiting  value  (Ml=.1907).  This  indicates  that  the  maximum 
opening  shock  force  occurs  during  the  over  inflation  phase  of 
deployment.  Calculation  of  the  opening  shock  force  at  t=tG 
produces  shock  magnitudes  in  excess  of  the  suspension  line 
constructed  strength  resulting  in  system  failure  prior  to  reaching 
the  calculated  tc*  The  parachute  of  example  1  is  not  a  practical 
application  for  the  heavier  weight  of  example  2. 
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SUMMARY  AND  CONCLUSIONS 

A  method  of  calculating  the  opening  shock  force  of  vertically 
deployed,  toward-the-earth  parachutes  has  been  developed.  The 
analysis  utilizes  a  dynamic  drag-area  signature  of  the  form  C[)S/Cd 
Sc  =  (1-r)  (t/t0)j  +  r  and  is  applicable  to  several  types  of 

parachutes.  Unlike  the  closed  form  analysis  of  Reference  1,  a 
furnished  computer  program  is  required  for  solution.  The  program 
has  two  modes  of  operation.  Mode  1  is  for  the  solid  cloth  canopy 
family  of  parachutes.  The  particular  inflation  reference  time  for 
the  system  and  deployment  condition  is  computed,  and  this  value  is 
used  to  calculate  the  instantaneous  drag  area,  velocity,  force 
profile,  and  distance  of  fall  during  the  inflation  process.  Mode 
2  requires  inflation  reference  time  as  input  data,  and  is  used  for 
all  types  of  parachutes.  This  mode  is  convenient  for  conducting 
surveys  to  evaluate  expected  opening  shock  force  variation  as  a 
function  of  inflation  time.  Vertical  opening  shock  forces 
calculated  by  this  technique  exceed  the  rule  of  thumb.  A  time 
step  of  dt=tQ/ 100000  in  the  program  gives  acceptable  results. 

The  inflation  reference  time  of  a  parachute  system  in 
vertical  fall  is  less  than  the  inflation  time  of  the  same  system 
deployed  horizontally  at  the  identical  altitude  and  velocity. 

The  inflation  reference  time  of  a  given  parachute  system,  for 
constant  launch  conditions,  depends  upon  the  system  weight.  For 
low  weight  finite  mass  assemblies  the  inflation  reference  time  is 
extended  due  to  the  trajectory  velocity  reduction  during 
inflation.  The  inflation  reference  time  reduces  rapidly  as  the 
system  weight  increases.  After  the  finite  mass  to  intermediate 
mass  transition  most  of  the  inflation  reference  time  decrease  has 
been  realized.  As  the  system  approaches  the  infinite  mass 
condition  the  inflation  reference  time  approaches  a  constant 
value . 


Parachutes  of  light  construction,  such  as  personnel  designs, 
deper  upon  the  reduction  of  trajectory  velocity  during  inflation 
for  i  ccessful  operation.  As  the  system  weight  is  increased  the 
Vt"  ry  reduction  is  not  realized  and  the  opening  shock  forces 
w_  eventually  exceed  the  parachute  constructed  strength, 
resulting  in  failure. 

A  closed  form  inflation  reference  time  equation  for 
horizontal  deployment  is  presented.  This  equation  demonstrates 
that  the  inflation  distance  (VstD)  is  a  function  of  altitude, 
canopy  cloth  permeability,  and  system  mass  and  area  ratios.  The 
particular  inflation  reference  time  solution  for  example  (1)  is  in 
close  agreement  with  the  empirical  approach.  While  more 
complicated,  the  theory  accounts  for  more  variables  and  is 
worthwhile  as  in  the  example  of  the  use  of  the  three  momme  silk 


j  > 
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cloth  at  sea  level  and  20/000  feet.  Usually,  the  inflation 
distance  is  considered  only  as  a  velocity-time  effect.  The 
closeness  of  the  theoretical  and  empirical  inflation  time  values 
at  low  altitudes  indicates  that  the  method  of  application  of  the 
following  parameters,  in  the  filling  time  calculation,  are 
reasonable: 

a)  Use  of  the  trajectory  velocity  as  the  canopy  inflow 
velocity 

b)  Use  of  the  cloth  permeability,  P,  to  determine 
instantaneous  cloth  airflow  rates 

c)  Correlation  of  the  instantaneous  mouth  inflow  area,  and 
instantaneous  pressurized  canopy  surface  area  with  the  drag  area 
ratio  variation 

d)  Canopy  volume  of  air,  1£0,  which  is  to  be  collected  in  the 
inflation  process. 

As  the  cloth  rate  of  airflow  exponent  "n"  increases  from  0.5, 
the  rate  of  cloth  airflow  increases  and  the  canopy  inflation  time 
is  extended. 
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Ab  3 1  rac  t 

An  analytical  method  of  calculating 
parachute  opening-shock  forces  based  upon 
wind-tunnel  derived  drag  area  time  signa¬ 
tures  of  several  solid  cloth  parachute 
types  i  r.  conjunction  with  a  scale  factor 
and  retardation  system  steady-state  param¬ 
eters  has  been  developed.  Methods  of 
analysing  the  inflation  time,  geometry, 
cloth  airflow  properties  and  materials 
elasticity  are  included.  The  effects  of 
mass  ratio  and  altitude  on  the  magnitude 
and  time  of  occurrence  of  the  maximum 
opening  shock  are  consistent  with  observed 
field  test  phenomena. 

I.  Introduction 

In  196*5 1  the  '(aval  Ordnance  Laboratory 
(riOL)  was  engaged  in  a  project  which 
utilized  a  35-foot-diameter ,  10-percent 
extended-skirt  parachute  (type  T-10)  as 
the  second  stage  of  a  retardation  system 
for  a  250-pound  payload.  Deployment  of 
the  T-10  parachute  was  to  be  accomplished 
at  an  altitude  of  100,000  feet.  In.  this 
rarefied  atmosphere,  the  problem  was  to 
determine  the  second  stage  deployment 
conditions  for  successful  operation,.  A 
search  of  available  field  tost  information 
indicated  a  lack  of  data  on  the  use  of 
solid  cloth  parachutes  at  altitudes  above 
30,000  feet. 

The  approach  to  this  problem  was  as 
follows:  Utilizing  existing  wind-tunnel 
data,  low-altitude  field  test  data,  ar.d 
reasonable  assumptions,  a  unique  engi¬ 
neering  approach  to  the  inflation  time 
and  opening-shock  problem  was  evolved  that 
provided  satisfactory  results.  Basically, 
the  method  combines  a  wind-tur.nel  derived 
drag  area  ratio  signature  as  a  function 
of  deployment  time  with  a  scale  factor  ar.d 
Newton's  second  lav;  of  motion  to  analyze 
the  velocity  and  force  profiles  during 
deployment.  The  parachute  deployment 
sequence  is  divided  into  tv/o  phases.  The 
first  phase,  called  "unfolding  phase," 
where  the  car.opy  is  undergoing  changes  in 
shape,  is  considered  to  be  inelastic  as 
ih°  pare.chut*-  \r  fle*es  initially  to  its 
steady-state  aerodynamic  size  for  the 
first  time.  At  this  point,  the  "elastic 
phase"  is  entered  where  it  is  considered 
that  the  elasticity  of  the  parachute 
materials  enter.-,  the  problem  -and  resists 
applied  forces  until  the  canopy  has 
r-'iehe-'i  full  Inflation. 

Th‘--  developed  equations  are  ir.  agree - 
r  t  v; :  th  the  observed  performance  of 
:;c. lid  ln*h  pari  *,vt°s  in  the  field,  r. u ■" :  1 
as  *.ue  f|.  ■r-eas'-  of  inflation  *:n-  or. 


altitude  increases,  effects  of  altitude  c>n 
open  '  rg-shock  force,  f  l:  it  ••  a.-d  i-finll-; 
mass  oparatio",  ar  d  inflat  io  dist.mce. 

II.  Development  of  Velocity  Patio  -and 
r'crce  Ratio  Equations  During 
the  Unfolding  Phase  of 
_ Parachute  Deployment _ 

The  parachute  deployment  v.ould  take 
clace  ir.  a  horizontal  attitude  in  accord¬ 
ance  with  Newton's  second  law  of  motion. 

£F  =  ma 


1 

2 


p'v?cds 


W  dV 
g  dt 


It  was  recognized  that  other  factors,  such 
as  included  air  mass,  apparent  mass,  and 
their  derivatives,  also  contribute  forces 
acti-  1  or.  the  system.  Since  definition 
of  t..ese  parameters  was  difficult,  the 
analysis  was  conducted  in  the  simplified 
form  shown  above.  Comparison,  of  calculated 
results  and  test  results  indicated  that 
the  omitted  terms  have  a  small  effect . 


Multiplying  the  right-hand  side  of  equa¬ 
tion  '(1)  by 


and  rearranging 


Ir,  ord^r  to  integrate  the  left-hand  term 
of  equation  (?),  the  drag  ar -a  ratio  must, 
re  defined  for  the  typ*>  of'  parachute  under 
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analysis  as  a  function  of  deployment 
i'j*Vyc:ico  time,  t0. 


FIC.  1  TYPICAL  INFINITE  MASS  FORCE-TIME  HISTORY  OF  A 
•  SOLID  CLOTH  PARACHUTE  IN  A  WIND  TUNNEL 
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REPRODUCED  FROM  REFERENCE  (II 


FIG.  3  TYPICAL  FORCE-TIME  CURVE  FOR  A  10%  EXTENDED 

SKIRT  PARACHUTE  UNDER  INFINITE  MASS  CONDITIONS. 
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REPRODUCED  FROM  REFERENCE  ID 
FIG.  4  TYPICAL  FORCE-TIME  CURVE  FOR  A  PERSONNEL 
GUIDE  SURFACE  PARACHUTE  UNDER  INFINITE 
MASS  CONDITIONS 


Figure  1  illustrates  a  typical  solid 
cloth  parachute  wind-tunnel  infinite  mass 
force- time  history  after  snatch.  In 
infinite  mass  deployment,  the  maximum  size 
ar.d  maximum  shock  force  occur  at  the  time 
of  full  inflation,  tf .  Hov.-ever,  tf  is 
inappropriate  for  analysis  since  it  is 
dependent  upon  the  applied  load,  structural 
strength,  and  materials  elasticity.  The 
reference  time,  tQ,  where  the  parachute  has 
attained  its  steady-state  aerodynamic  size 
for  the  first  time,  is  used  as  the  basis  for 
performance  calculations. 

At  any  instant  during  the  unfolding 
phase,  the  force  ratio  F/Fs  can  be  deter¬ 
mined  as  a  function  of  the  time  ratio, 
t/tc* 


-  p  V  cC-,S 
~  s  Do 


Since  the  wind-tunnel  velocity  and  density 
are  constant  during  infinite  mass  deployment 


CDS 

CDSo 
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FIG.  2  TYPICAL  FORCE-TIME  CURVE  FOR  A  SOUO  FLAT  PARA¬ 
CHUTE  UNOER  INFINITE  MASS  CONDITIONS. 


FIG.  5  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  ELLIPTICAL 
PARACHUTE  UNOER  INFINITE  MASS  CONDITIONS 


FIG.  6  TYPICAL  FORCE-TIME  SIGNATURE  FOR  THE  RING 

SLOT  PARACHUTE  20%  GEOMETRIC  POROSITY  UNOER 
INFINITE  MASS  CONDITIONS 


Infinite  mass  opening- shock  signatures 
cf  several  types  of  parachutes  arc  pre¬ 
sented  in  Figures  2  through  6.  Analysis 
of  these  signatures  using  the  force  ratio, 
7/?.,  -  time  re.tio,  t/t0,  technique 
indicated  a  similarity  in  the  performance 
of  the  various  solid  cloth  types  of 
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parachutes  which  wore  examined.  The  geo¬ 
metrically  porous  ring  slot  panchute 
displayed  a  completely  different  signature, 
as  was  expected.  These  data  are  illus¬ 
trated  in  Figure  7.  If  an  initial  boundary 


FIG.  7  DRAG  AREA  RATIO  VS.  TIME  RATIO 


condition  of  CDS/CDSc  =  0  at  time  t/tQ  =  0 

is  assumed,  then,  the  data  car.  be  approxi¬ 
mated  by  fitting  a  curve  of  the  form 


At  the  time  that  equation  ('•)  was  ascer¬ 
tained,  it  suggested  that  the  geometry  of 
the  deploying  parachute  independent 

of  density  and’ velocity.  It  was  also 
postulated  that  although  this  expression 
had  been  determined  for  the  infinite  mass 
condition,  it  would  also  be  true  for 
the  finite  mass  case,  '.his  phenomenon 
has  since  been  independently  observed  and 
confirmed  by,  Berndt  and  lie  ’..‘eese  in 
reference  ($). 

Since  the  drag  area  ratio  was  determined 
from  actual  parachute  deployments,  it  was 
assumed  that* the  effects  of  apparent  mass 
and  included  mass  or.  the  deployment  force 
history  were  accommodated. 

The  right-hand  term  of  equation  (2) 
contains  the  expression 


*SVoCDSo 


This  term  can  be  visualised  as  shown  in 
Figure  8  to  be  a  ratio  cf  the  retarded 
mass  (including  the  parachute)  to  an 
associated  mass  of  atmosphere  contained 
in  a  right  circular  cylinder  -which  is 
generated  by  moving  an  inflated  parachute 
of  area  CjSq  for  a  distance  equal  to  the 
product  of  7st0  through  an  atmosphere  of 
density,  o. 


A  more  realistic  drag  area  ratio  expres¬ 
sion  was  determined  which  includes  the 
effect  of  initial  area  at  line  stretch. 


CDS 

CDSo 


FIG.  8  VISUALIZATION  OF  THE  MASS  RATIO  CONCEPT 


The  mass  ratio,  M,  is  the  scale  factor 
which  controls  the  velocity  and  force 
profiles  during  parachute  deployment. 
Substituting  it  and  into  equa¬ 

tion  (2),  integrating,  and  solving  for 

y/v. 


where  n  is  the  ratio  of  the  projected 
mouth  area  at  line  stretch  to  the  nteidy- 
ctftte  projected  frontal  area.  Expandin', 
equation  (4) 
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The  instantaneous  shoe);  factor  is  defined 

as 

f 

X.  »  —  -  * - - - 

T*  |pV82cdSo 


Figures  9  and  10  illustrate  the  velocity 
and  force  prof  lice  .generated  fror.  cquatUnn 
(7)  and  (8)  for  initial  projected  area 
ratios  of  n  -  0,  and  0.!?  with  various  mass 
ratios. 

IV.  Methods  for  Calculation  of 
the  Reference  Time,  t0 


If  the  altitude  variation  during  deployment 
is  snail,  then,  the  density  nay  be  consid¬ 
ered  os  constant 


from  equations  (5)  and  (7) 


The  ratio  concept  ic  an  ideal  method  to 
analyze  the  effects  of  the  various  param¬ 
eters  on  the  velocity  and  force  profiler, 
of  the  opening  parachutes;  however,  a 
means  of  calculating  tD  is  required  to  fore 
specific  values  can  be  computed.  Methods 
for  computing  the  varying  mass  flcv;  i  r.to 
the  inflating  canopy  mouth,  the  varying 
mass  flow  out  through  the  varying  inflated 
canopy  surface  area,  and  the  volume  c-f 
air,  which  must  be  collected  during 
the  inflation  process  are  required. 


III.  Maximum  Shock  Force  and  Time  of 
Occurrence  During  the  Unfolding  Phase 

The  time  of  occurrence  of  the  maximum 
Instantaneous  shock  factor,  x*,  is  diffi¬ 
cult  to  determine  for  the  general  ease. 
However,  for  r  =  0,  the  maximum  shock 
factor  and  time  of  occurrence  are  readily 
calculated.  For  n  **  0 


Figure  11  represents  a  solid  cloth- 
type  parachute  canopy  at  some  instant 
during  inflation.  At  any  given  instant, 
the  parachute  drag  area  is  proportional 
to  the  maximum  inflated  diameter.  Also, 
the  maximum  diameter  in  conjunction  with 
the  suspension  lines  determines  the  Inflow 
mouth  area  (A-A)  and  the  pressurized 
canopy  area  (B-B-B).  This  observation 
provided  the  basis  for  the  following 
assumptions.  The  actual  canopy  shape  is 
cf  minor  importance. 

a.  The  ratio  of  the  instantaneous 
mouth  inlet  area  to  the  steady-state 
mouth  area  is  in  the  same  ratio  as  the 
instantaneous  drag  are-a. 


Setting  the  derivative  of 
to  time  equal  to  zero  and 
t/tQ  at  x*  m&x 


with  respect 
solving  for 


(9) 


and  the  maximum  shock  factor  is 


6 

„  16  /21M\7 
1  *"  49  \  4  ) 


(10) 


_  CDS 
^o  CDSo 

b.  The  ratio  of  the  instantaneous 
pressurized  cloth  surface  area  to  the 
canopy  surface  area  is  in  the  same  ratio 
as  the  instantaneous  drag  area. 


CDS 


So  Vo 


c.  Since  the  suspension  lines ir.  the 
unpressurized  area  of  the  canopy  are 
straight,  a  pressure  differential  has  not 
developed,  and,  therefore,  the  net  air¬ 
flow  in  this  zone  is  zero. 

Based  on  the  foregoing  assumptions, 
the  mass  flow  equation  can  be  written 


Equations  (9)  and  (10)  are  valid  for  value* 
of  K  *  ~~~  (0.19),  since  for  larger  values 

of  ti,  the  maximum  shock  force  occurs  in 
the  elastic  phase  of  inflation. 


dm  »■  m  inflow  -  m  outflow 


dv 

dt 


Av\: 


pa0p 


A-B 


>w  v*.  w  ia.  i  Jw-tnri/vuvuvWLW* 


V/V  VELOCITV  RATIO  V/V.  VELOCITY  RATIO 
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dv  c-  s  c.,i; 

p  —  =  pVA,.  — - pA,  — —  p  (ill 

H  1.0  ,,  ,  rt,0  \ 


"r°k  (^)  1(0  iTi  av 


o'  J  (1? 


I integrating: 


FIG.  11  PARTIALLY  INFLATED  PARACHUTE  CANOPY 


From  equation  (3) 


fey 


;  for  ri  =  0 


CDSC  Uo 


From  equation  (7) 


— 2 - ;  n  =  0 

1  /t  \? 


■'rom  equation  (Po) 


‘  (?r 


V  =  A,,  V  t  C.  tr, 

—  O  !to  S  O 


hi] 


/c^\nl/t\c  f  v-  Y 

*s°kU/  (J  1  +  i/r\7 

Jo  1  VM  Vt0/  J 


Measured  values  of  n  Indicate  a  data 
range  from  0 through  0.771.  A 
convenient  solution  to  the  reference  time 
equation  evolves  when  n  Is  assigned  a 
value  of  1/?.  Integrating  equation  (13) 
and  using 


vat  (•:  = 

P  O 


«pCD;:o 


LET  Y*x  = 


V  -  Aso*  f)‘ 


«WdSo 


Equation  (14)  expresses  the  unfolding 
reference  tine,  t0,  in  terms  of  mass, 
altitude,  snatch  velocity,  airflow  char¬ 
acteristics  of  the  cloth,  and  the  steady- 
state  parachute  geometry.  Note  that  the 
term  gp'.'r/V;  is  the  ratio  of  the  included 
air  mass  to  the  mass  of  the  retarded 
hardware.  Multiplying  both  sides  of 
equation  (14)  by  Ve  demonstrates  that 


^  “  Ws 


1  +  -- 


- ,  d  t 


\„t  -  a  constant  which  Is 

°  0  a  fur.ctlor,  of  altitude 


Figures  IP  and  lj  indicate  tic-  para¬ 
chute  unfolding  time  and  Unfolding 
distance  for  valu'-s  of  n  -  1/;  and  n  - 
O.A3pl:f,.  Note  the  variation  and  co.ver- 
g'Tice  1 1 •  i  rising  altitude.  opeiilr.,-- 

siioeF  fore--  is  strongly  ii.flu'  n.  ■'  ;  by 
the  infJ  at  lot:  tlr..-.  Neoause  of  ‘  1,  i  n. ,  th' 


■  jrA'WAir.P.»PWX  IFm  VPi.’ITKIt ' 
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\ 


! 


value  or  tc  calculated  by  using  a.  real¬ 
istic  value  of  r.  should  be  used  in  the 
lower  atmosphere . 

As  an  example  of  this  method  of  opening- 
shock  analysis,  let  us  examine  the  effect 
of  altitude  on  the  openi hr- shock,  force  of 
a  T-10-type  parachute  retarding  a  P00- 
pound  weight  from  a  snatch  velocity  of 
V8  =  400  feet  per  second  at  sea  level. 
Conditions  of  constant  velocity  ar.d  con¬ 
stant  dynamic  pressure  are  investigated. 

The  results  are  presented  in  Figure  14. 

At  low  altitudes,  the  opening-shock  force 
is  lesc  than  the  steady-state  drag  force: 
however,  as  altitude  rises,  the  opening 
shock  eventually  exceeds  the  steady-state 
drag  force  at  some  altitude.  This  trend 
is  in  agreement  with  field  test  observations. 

V.  Correction  of  tc  for 
Initial  Area  Effects 


The  unroldlng  reference  time,  tD, 
calculated  ,,y  the  previous  methods  assumes 
that  the  parachute  inflates  from,  zero  drag 
area.  In  reality,  a  parachute  has  a  drag 


TIG  t7  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  TIME 
"I  ■■  AT  CONSTANT  VELOCITY  AND  CONSTANT 
DYNAMIC  PHESSURE  FOR  n  •  1/7  ANO  n  ■  063796 


CLOTH  MIL  C  7070.  TYPE  III 

36  FOOT.  1C  .  EXTENDED  SKIRT  CANOPY 


HG  >3  EFFECT  OF  ALTITUDE  ON  THE  UNFOLDING  OIS 

TANCE  AT  CONSTANT  VELOCITY  AND  CONSTAN1 
DYNAMIC  PRESSURE  FOR  n  »  1/7  AND  n  •  0.63246. 


FIG  14  VARIATION  Of  STEADY  STATE  DRAG.  Ff.  ANO  MAXI 

MUM  OPENING  SHOCK  WITH  AITITUPE  FOR  CONSTANT 
VELOCITY  AND  CONSTANT  DYNAMIC  PRESSURE 


area  at  the 
•.r,  has  beer, 
'  applied, 

the  initial 


begin'- ing  of  inflation, 
calculated,  a  correction 
based  upon  what  Is  known 
cor  d  1 t  ions . 


Once 

C-ar 

aoout 


AH 
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Case  A  -  '.inn 
; re  a  is  known 


t  r. :  *  Lai  pro  Joe  tec 


than  the  maximum  extensibility,  c  ,  o:' 
the  materials .  '!llX 


1/3 


calculated 


Development  of  the  Analysis  1'  -he 
elastic  phase  of  inflation  is  sir'lar  to 
the  technique  used  ir  the  unfoldirg  phase 
Hewton's  second  la1.;  of  .to  f  L  or.  is  us-d, 
together  with  the  drag  area  ratio  signa¬ 
ture  and  rr.ass  ratio 


CDSo 


corrected 


,V3' 


calculated 


(iz) 


which  Is  still  valid, as  show a  in  Figure  7 

t  V 

/  »  \6 

dt  =  V 


Case  B  -  Vine:  the  initial  drag  area 
Is  known 


'D  i 


CDSo 


/  c  s  \ 

"  ■(&) 


calculated 


Integrating  and  solving  for  — 

Vs 


V0  7M 


(t)7- 


(17) 


v.nere  —  is  the  velocity  ratio  of  the 


corrected 


1/61 


ur. folding  process  at  tir.e  t  =  tc- 


ealeu lated 

(16) 


Trie  mass  ratio  should  'cw  be  adjusted  for 
the  correced  t0  before  velocity  a-d  force- 
profiles  arc  deterr.i- ed. 

VI.  Open  1 ' g-3hcck  Force,  Velocity  Ratio, 
and  I'flation  Tine  During  the 
Elastic  Phase  of  parachute  Inflation 

me  mass  ratio,  !■',  is  an  important 
parameter  ir  parachute  analysis.  For 
values  of  M  <<4/21,  the  maximum  opening- 
shock  force  occurs  early  ir.  the  Inflation, 
process,  and  the  elastic  properties  of 
the  canopy  are  rot  significant.  As  the 
rass  ratio  approaches  =  4/51,  the 
magnitude  of’  the  cp-r  Irg-shock  force 
increases,  ad  the  time  of  occurrence 
happens  later  in  the  deployment  sequence. 
For  mass  ratios  >  4/21,  the  maximum 
shock  force  will  occur  after  the  reference 
*  i  ir.e ,  tc.  Parachutes  designed  for  high 
rass  ratio  c-peration  must  provide  a 
structure  cf  sufficient  co-structed 
strength,  Fc ,  sc  that  the  actual  elon¬ 
gation  of  the  car-cpy  under  load  is  less 


V 

o  _ _ 

1 

Vs  i  +  i 

11  -  !l1?  +  rll-r.  "il  +  ^1 

M 

7  ?  J 

( 

.he  instantaneous  shock  factor  In  the 

elastic  prase 

becomes 

CDSo  ' ’ s 


H.  (•)'• 


(t/i 


The  end  point  of  the  inflation  process 
depends  upon,  the  applied  loads,  elasticity 
of’  the  ca-cpy,  and  the  cor.  struct  el  strength, 
of  the  car-uhute.  A  linear  load  ele  cation 
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relationship  is  utilized  to  determine  the  The  next  force  ir.  the  series  at  constant  q 
maximum  drag  area. 

-  v  A] 

c  £r.ax  r  1  -  ors  7  " 

—  =  — — — 


=  0  -  Cn)‘ 


The  force,  F,  is  initially  the  instanta¬ 
neous  force  at  the  end  of  the  unfolding 
process 


F  =  X  F 
o  s 


Subsequent  elongations  lr.  the  system  can 
be  shown  to  be 


C1  =  co  O  J  €o)‘ 


where  X0  is  the  shock  factor  of  the 
unfolding  phase  at  t  =  tQ 


*o  =  r — ; — 

i  +  i 

H  l  7 


-_nj_  +  r!(l_-_n_)  +  ^2 
7  2 


Since  the  inflated  shape  Is  defined,  the 
drag  coefficient  is  considered  to  be 
constant,  and  the  instantaneous  force  is 
proportional  to  the  dynamic  pressure  and 
projected  area.  The  maximum  projected 
area  would  be  developed  if  the  dynamic 
pressure  remained  constant  during  the 
elastic  phase.  Under  very  high  mass 
ratios,  this  ia  nearly  the  case  over  this 
very  brief  time  period;  but  as  the  mass 
ratio  decreases,  the  velocity  decay  has 
a  more  significant  effect.  The  simplest 
approach  for  all  mass  ratios  is  to 
determine  the  maximum  drag  area  of  the 
canopy  as  if  elastic  inflation  had 
occurred  at  constant  dynamic  pressure. 

Then  utilizing  the  time  ratio  determined 
as  an  end  point,  intermediate  shock 
factors  can  be  calculated  from  equation 
(19)  and  maximum  force  assessed. 

The  initial  force,  X0FS,  causes  the 
canopy  to  Increase  in  projected  area.  The 
new  projected  area  in  turn  Increases  the 
total  force  on  the  canopy  which  produces 
a  secondary  projected  area  increase.  Tr.e 
resulting  series  of  events  are  resisted 
by  the  parachute  materials.  Tr.e  parachute 
must,  therefore,  be  constructed  of  suffi¬ 
cient  strength,  to  prevent  the  elongation 
of  the  materials  from  exceeding  the 
maximum  elongation. 


XoFs 

c  - - e 

o  P  max 


£?  =  co  (!  +  £o  (X  +  £o)  )' 

The  required  canopy  constructed  strength 
can  be  determined  for  a  giver  set  of 
deployment  conditions.  The  limiting  value 
of  the  series  (*t)  determines  the  end 
point  time  ratio. 

(l£\6  a  =  (1  +  c  ){? 

Uo  /  CDSo 


(t)W:  = 


(i  +  t,y 


Figure  15  illustrates  the  maximum  drag 


area  ratio  as  a  function  of  e0. 
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F'G.  15  MAXIMUM  DRAG  AREA  RATIO  VS  INITIAL 
SLONCATION 
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VII.  Application  of  Cloth  Permeability 

to  the  Calculation  of  the  Inflation 
Time  of  Solid  Cloth  Parachutes _ 

The  mass  outflow  through  the  presBur- 
»r  %n  lnra«tlnR  •alia  cloth 
parachute  at  any  instant  is  dependent  upon 
the  canopy  area  which  is  subjected  to  air¬ 
flow  and  the  rate  of  airflow  through  that 
area.  The  variation  of  pressurized  area  os 
a  function  of  reference  time,  te,  was 
earlier  assumed  to  be  proportional  to  the 
Instantaneous  drag  area  ratio,  leaving  the 
rate-of- airflow  problem  to  solve.  The 
permeability  parameter  of  cloth  was  a 
natural  choice  for  determining  the  rate  of 
airflow  through  the  cloth  as  a  function  of 
pressure  differential  across  the  cloth. 
Heretofore,  these  data  have  been  more  of 
a  qualitative,  rather  than  quantitative, 
value..  A  new  method  of  analysis  was 
developed  wherein  a  generalized  curve  of 
the  form  P  =>  k(AP)n  was  fitted  to  cloth 
permeability  data  for  a  number  of  different 
cloths  and  gives  surprisingly  good  agree¬ 
ment  over  the  pressure  differential  range 
of  available  data.  The  pressure  differ¬ 
ential  was  then  related  to  the  trajectory 
conditions  to  give  a  generalized  expression 
which  can  be  used  in  the  finite  mass  ratio 
range,  as  well  as  the  infinite  mass  case. 
The  permeability  properties  were  trans¬ 
formed  into  a  mass  flow  ratio,  M',  which 
shows  agreement  with  the  effective  porosity 
concept. 


Measured  and  calculated  permeability 
pressure  data  for  several  standard  cloths 
are  illustrated  in  Figure  16.  This  method 
has  been  applied  to  various  types  of  cloth 
between  the  extremes  of  a  highly  permeable 
3-momme  silk  to  a  relatively  impervious 
parachute  pack  container  cloth  with 
reasonably  good  results,  see  Figure  17. 


The  car.opy  pressure  coefficient,  Cp, 
is  defined  as  the  ratio  of  the  pressure 
differential  across  the  cloth  to  the 
dynamic  pressure  of  the  free  stream. 


C  _  AP  _  P( internal)  -  P(external) 
p  q  1/2  pV2 

(25) 

where  V  is  based  on  equation  (7). 

The  permeability  expression,  P  =>  k(dP)n 
becomes 

2  n 

?  -  -(Cp  (25) 


Although  some  progress  has  been  made 
by  Melzig  and  others  on  the  measurement  of 
the  variation  of  the  pressure  coefficient 
on  an  actual  inflating  canopy,  this 
dimension  ar.d  its  variation  with  time  are 
still  dark  areas  at  the  time  of  this 
writing.  At  the  present  time,  a  constant 
average  value  of  pre' jure  coefficient  is 


UBed  in  these  calculations .  Figure  ifi 
presents  the  effect  ot  pressure  coefficient 
ar.d  altitude  on  the  unfolding  time  for 
constant  deployment  conditions. 

It  le  wall  Known  that  the  Inflation 
time  of  solid  cloth  parachutes  decreases 
as  the  operational  altitude  increases. 

This  effect  can  be  explained  by  consid¬ 
ering  the  ratio  of  the  mass  outflow 
through  a  unit  cloth  area  to  the  mass 
inflow  through  a  unit  mouth  area. 


M'  =s  mass  flow  ratio 


mass  outflovr 
maos  inflow 


where 


mass  outflow  -  P-  — (per.  ft2 cloth  o.rea) 
r  ft  -sec 


and 


mass  inflow  =  V  -~HS5-  (per  ft2  inflow  area) 

P  ft^-Bec 


REPRODUCED  FROM  REFERENCE  (4) 

FID.  16  NOMINAL  POROSITY  OF  PARACHUTE  MATERIAL 
VS  DIFFERENTIAL  PRESSURE. 
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Ircrefore,  the  r.ass  flow  ratio  becomes 


Effect lve  porosity,  C,  Is  defined  as 
the  ratio  of  the  velocity  through  the 
cloth,  u,  to  a  fictitious  theoretical 
velocity,  v,  v.'hic'^  will  produce  the  par¬ 
ticular  ap=  i/?ov  ■ . 


fIC  17  COMPARISON  Of  MEASURED  ANO  CALCULATED 

PERMEABILITY  FOR  RELATIVELY  PERMEABLE  ANO 
IMPERMEA01  E  CLOTHS 


AVERAGE  CANOPY  PRESSURE  COEFFICIENT 
DURING  INFLATION  INCLUDING  THE  VENT 


FIG  IB  EFFECT  OF  PRESSURE  COEFFICIENT  AND  ALTITUDE 
ON  THE  UNFOLDING  TIME 


effective  porosity,  C  =  —  (T3) 


Comparison  of  the  mass  flc1.'  ratio  and 
previously  published  effective  porosity 
data  is  shown  lr,  Figure  19.  The  effects 
of  altitude  ar.d  velocity  or.  the  r.ass  flov; 
ratio  are  presented  ir.  Figures  20,  ar.d  21 
for  constant  .elocity  ar.d  constant 
alt'tude.  Tits  decrease  c.f  cloth  perme¬ 
ability  •..•1th  altitude  Is  evident. 

The  permeability  cor.sta- ts  1  V  and  "r." 
can  be  determined  frorr.  the  per.teabl  11  ty 
pressure  differential  data  as  obtained 
frc~  ar.  ir  strur.e:  t  such  as  a  Frasier 
.-crr.ear.eter.  Two  data  points,  "A1'  and 


REPRODUCED  FROM  REFERENCE  Ml 


TIG  19  THE  EFFECTIVE  PO  IOSITY  OF  PARACHUTE 
MATERIALS  VS.  DIFFERENTIAL  PRESSURE 
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"B,  "  are  selected  ir.  such  a  manner  that 
point  "A"  is  in  a  low-pressure  zone  below 
the  knee  of  the  curve,  and  point  3  is 
located  in  the  upper  end  of  the  high- 
pressuro  zone,  as  shown  in  Figure  22. 

The  two  standard  measurements  of  1/2 
inch  of  v.-ater  and  20  inches  of  water 
appear  to  be  good  data  points  if  both  are 


FIG.  20  EFFECT  OF  ALTITUOE  ON  MASS  FLOW  RATIO  AT 
CONSTANT  VELOCITY 


FIG.  21  EFFECT  OF  VELOCITY  ON  MASS  FLOW  RATIO  AT 
CONSTANT  DENSITY 


FIG.  22  LOCATION  OF  DATA  POINTS  FOR  DETERMINATION 
OF  ”k“  AND  “n" 


available  on  the  same  sample.  Substi¬ 
tuting  the  data  from  points  A.  and  "B 

into  P  »  k(AP)n  : 


VIII.  Determination  of  the 
Parachute  Included  Volume 
and  Associated  Air  Mass 


Before  the  reference  time,  tc,  and 
inflation  time,  tf,  can  be  calculated, 
the  volume  of  atmosphere,  Vo ,  which  is 
to  be  collected  during  the  inflation 
process  must  be  accurately  known.  This 
requirement  dictates  that  a  realistic 
Inflated  canopy  shape  and  associated 
volume  of  atmosphere  be  determined. 

Figure  23  was  reproduced  from  reference 
(5).  The  technique  of  using  lampblack 
coated  plates  to  determine  the  airflow 
patterns  around  metal  models  of  inflated 
canopy  shapes  was  used  by  the  investigator 
of  reference  (5)  to  study  the  stability 
characteristics  of  contemporary  parachutes, 
i.e.,  1943.  A  by-product  of  this  study 
is  that  it  is  clearly  shown  that  the 
volume  of  air  within  the  canopy  bulges 
out  of  the  canopy  mouth  (indicated  by 
arrov;s)  and  extends  ahead  of  the  canopy 
•_  hem.  This  volume  must  be  collected  during 
the  inflation  process.  Another  neglected, 
but  significant,  source  of  canopy  volume 
exists  ir.  the  billowed  portion  of  the  gore 
panels . 
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VENT PARACHUTE 


REPRODUCED  FROM  REFERENCE  (SI 

FIG.  23  AIRFLOW  PATTERNS  SHOWING  AIR  VOLUME 
AHEAD  OF  CANOPY  HEM 


The  steady-state  canopy  shape  has  been 
observed  in  wind-tunnel  and  field  tests  to 
be  elliptical  in  profile.  Studies  of  the 
inflated  shape  and  included  volume  of 
several  parachute  types  (flat  circular, 

10  percent  extended  skirt,  elliptical, 
hemispherical,  ring  slot,  ribbon,  and 
cross)  are  documented  in  references  (6) 
and  (7).  These  studies  demonstrated  that 
the  steady-state  profile  shape  of  inflated 
canopies  of  the  various  types  can  be 
approximated  to  be  two  ellipses  of  common 
major  diameter,  2a,  and  dissimilar  minor 
diameters,  b  and  b' ,  as  shown  in  Figure  24. 
It  was  also  shown  that  the  volume  of  the 
ellipsoid  of  revolution  formed  by  revolving 
the  profile  shape  about  the  canopy  axis 
was  a  good  approximation  of  the  volume  of 
atmosphere  to  be  collected  during  canopy 
inflation  and  included  the  air  volume 
extended  ahead  of  the  parachute  skirt  hem 
together  with  the  billowed  gore  volume. 


FIG.  24  PARACHUTE  CROSS  SECTION  NOMENCLATURE 


Tables  I  and  II  are  summaries  of  test 

results  reproduced  from  references  (6)  and 

(7),  respectively,  for  the  convenience  of 

the  reader. 
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TABLE  I  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  12-GORE  AMD  16-GORE  CONFIGURATIONS 


Type 

- 1 

No  o* 

Gom 

Suip«ruion 

Or*  length 

mfhti 

Velocity 

n\ph  Ipi 

2a 

O 

o 

Scale  factor  K 
2i  2i 

r~iro  1 

N 

4 

A.t>  Rjlio 

b  b 

i  i 

b  b 

4  4 

Volume  tn 

Vm  Vc 

V 

V 

cr 

Flat  Circular 

12 

34 

SO 

73 

646 

660 

856 

6115 

8817 

1  4932 

4476 

4481 

6980 

1  66 

16 

34 

SO 

73 

663 

669 

820 

5568 

9039 

1  4597 

4450 

4100 

7325 

1  65 

10%  £mt*nd*d 

12 

34 

100 

147 

663 

65  2 

881 

6424 

8860 

1  5284 

3926 

4400 

6783 

1  73 

Skirl 

16 

34 

17 

25 

664 

640 

785 

5580 

850? 

1  4082 

4061 

3920 

6197 

1  53 

€  Hipt'c  ji 

12 

34 

7S 

110 

916 

812 

5626 

9657 

1  6283 

3322 

6406 

16 

34 

17 

26 

875 

800 

6169 

8163 

1  4332 

2726 

4405 

M*mi»ph«fK#l 

12 

34 

ITS 

183 

996 

1  254 

1  0005 

9080 

1  9085 

6224 

8606 

16 

34 

75 

1  10 

994 

1  186 

9129 

9380 

1  8509 

5921 

83  70 

ftmyjlot 

12 

34 

TS 

3? 

607 

664 

853 

6566 

8735 

i  630 

3800 

3660 

5903 

1  66 

16%  Gaomatf.C 

12 

34 

100 

147 

616 

663 

922 

6566 

8736 

1  530 

3800 

4198 

6166 

1  62 

Poronty 

12 

34 

200 

293 

637 

686 

918 

6566 

8736 

1  530 

3800 

4624 

6826 

1  90 

16 

34 

23 

37 

611 

658 

827 

6004 

8890 

1  4894 

3800 

3763 

5686 

1  60 

16 

34 

100 

147 

617 

664 

864 

6004 

8890 

1  4894 

3800 

3986 

6030 

1  59 

16 

34 

200 

293 

64S 

696 

844 

6004 

8890 

1  4894 

3800 

(- 

4430 

6897 

1  8  2 

Ribbon 

12 

34 

25 

37 

586 

632 

859 

6558 

8768 

1  5326 

3800 

3323 

5336 

1  40 

?4%  G«o<nrii.< 

12 

34 

100 

147 

615 

663 

837 

6568 

0768 

1  5326 

3800 

3714 

6163 

1  62 

Pofo'.'ty 

12 

34 

200 

293 

$32 

681 

877 

6658 

8768 

1  5326 

3800 

4280 

6683 

1  76 

16 

34 

25 

37 

603 

650 

797 

5570 

85  78 

1  4148 

3800 

3438 

5368 

1  41 

16 

M 

100 

147 

C26 

674 

791 

5570 

8578 

1  4148 

3800 

3804 

6983 

1  67 

16 

34 

200 

293 

648 

698 

781 

5570 

8578 

1  4148 

3800 

4164 

6656 

1  75 

C'OII  Chull 

34 

25 

37 

710 

543 

1  24* 

886  7 

1  27  76 

2  1643 

1928 

3768 

5798 

301 

*■'1  *  2*4 

34 

100 

147 

707 

540 

1  270 

0867 

1  2776 

2  1643 

1928 

3810 

6712 

2  96 

34 

2  00 

293 

716 

64  7 

1  285 

8867 

1  2776 

2  1643 

1928 

4  ?  1 2 

5976 

307 

47 

26 

37 

769 

680 

1  113 

8494 

1  261? 

2  1006 

1928 

405? 

6868 

3  56 

47 

100 

147 

>M 

65  7 

1  206 

8494 

1  2512 

2  1006 

1928 

3973 

5968 

3  09 

47 

t—..  .  .  .  ^ 

200 

293 

776 

69? 

1  110 

8494 

1  261? 

2  1006 

1928 

4292 

7303 

3  79 
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TABLE  II  SUMMARY  OF  PARACHUTE  SHAPE  TEST  RESULTS 
FOR  24-GORE  AND  30-GORE  CONFIGURATIONS 


P*»*chui» 

Typ* 

No  of 

Gom 

Sutptntion 

Line  Length 

■flfhri 

Veloc 

ly 

Ipi 

Scale  factor  K 

22  2i 

o'.  0, 

N 

4 

A 

b 

a 

«ei  Ratio 

r> 

a 

h  b 

a  *  a 

VM 

Volume  in  ^ 

V 

Vf 

\Jn 

vH 

f  l*i  Cucuiin 

24 

34 

50 

73 

677 

679 

795 

5  758 

0126 

1  3884 

436? 

4695 

7273 

1  67 

30 

34 

1  7 

26 

668 

669 

827 

6214 

7806 

1  4020 

4342 

4626 

702  7 

1  6? 

10%  E  *  1  mdvd  * 

24 

34 

100 

147 

666 

640 

834 

5949 

87  71 

1  4  770 

4138 

4446 

6930 

1  67 

Skirl 

30 

34 

1  7 

25 

650 

633 

825 

6255 

796/ 

14  177 

417? 

4076 

6266 

1  60 

Ring  Slot  16% 

24 

34 

25 

37 

663 

665 

824 

5800 

9053 

1  4853 

3591 

3878 

6031 

1  68 

Geometrically 

24 

34 

100 

147 

680 

682 

819 

5800 

9053 

1  4053 

3691 

4079 

6610 

1  01 

Poroirt 

24 

34 

200 

293 

694 

696 

809 

5800 

9053 

1  4863 

3691 

4270 

6924 

1  93 

30 

34 

26 

37 

677 

678 

708 

5800 

9053 

1  4863 

358? 

3826 

64  04 

1  79 

30 

34 

.00 

147 

684 

685 

807 

5800 

9053 

1  4863 

368? 

4023 

6580 

1  84 

JO 

34 

200 

293 

698 

699 

800 

5800 

9053 

1  4053 

3682 

4  260 

701? 

„e 

H  bbon  24\ 

24 

34 

25 

37 

671 

673 

7  70 

5980 

8187 

1  4167 

3691 

3591 

5968 

1  66 

GfOmft'iCjfly 

24 

34 

100 

147 

6  76 

678 

813 

5980 

8187 

1  4167 

3591 

392  7 

609  7 

1  70 

Porout 

24 

34 

200 

293 

68  7 

689 

804 

5980 

e  ie  7 

1  4  167 

3691 

4061 

6389 

1  79 

30 

34 

25 

37 

655 

65  7 

78? 

6021 

8463 

1  4484 

368? 

1396 

5666 

1  50 

JO 

* 

100 

147 

669 

670 

784 

6021 

8463 

1  4434 

.168? 

362? 

602? 

1  68 

30 

X 

20 

293 

677 

679 

823 

6021 

8463 

1  4404 

358? 

400? 

6256 

1  75 

‘Smir  thit  parachute  ah 

breathing  dor- 

nq  f hr  fMf 

ir«rr* 

photograph!  «*#/# 

taken  at 

each  ipr»4 

T  Hr  da  la 

Airir  rrdur 

ed  *  run 

th»  plwtnq 

aph  rffbifh  moil 

'Oionibfy  appeared  lo  iep'*l*nt  ihe  ttal* 
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X.  List  of  Symbols 

-  Steady- state  projected  area  of  the 
inflated  parachute,  ft2 

-  Instantaneous  canopy  mouth  area,  ft~ 


Ako  -  St|ady-state  Inflated  r.outh  area, 


a 

2a 


max 


g 

k 

m 

M 


-  Acceleration,  ft/sec^ 

-  Maximum  Inflated  parachute  diameter 
of  gore  mains earn,  ft 

-  Minor  axis  of  the  ellipse  bounded 
by  the  major  axis  (2a)  and  the  ver.t 
of  the  canopy,  ft 

-  Minor  axis  of  the  ellipse  which 
includes  the  skirt  hen  of  the 
canopy,  ft 

-  Effective  porosity 

-  Parachute  coefficient  of  drag 

-  Parachute  pressure  coefficient, 
relates  internal  and  external 
pressure  jaP)  on  canopy  surface  to 
the  dynamic  pressure  of  the  free 
stream 

-  Nominal  diameter  of  the  aerodynamic 

decelerator  -  ft 

-  Instantaneous  force,  lbs 

-  Steady- state  drag  force  that  would 
be  produced  by  a  fully  open  para¬ 
chute  at  velocity  Vg,  lbs 

-  Constructed  strength  of  the  para¬ 
chute,  lbs 

Maximum  opening-shock  force,  lbs 

o 

-  Gravitational  acceleration,  ft/sec 

-  Permeability  constant  of-  canopy 
cloth 

-  Mass,  slugs 

-  Mass  ratio  -  ratio  of  the  mass  of 
the  retarded  hardware  (including 
parachute)  to  a  mass  of  atmosphere 
contained  in  a  right  circular 
cylinder  of  length  (V6t0),  face  area 

(cdso)'  and  denEity  (o) 

-  Mass  flow  ratio  -  ratio  of  atmo¬ 
sphere  flowing  through  a  unit  cloth 
area  to  the  atmosphere  flowing 
through  a  unit  inlet  area  at 
arbitrary  pressure 

-  Permeability  constant  of  canopy 
cloth 


S  =A 
o  so 


Cloth  permeability  -  rat?  of  air¬ 
flow  through  a  cloth  at  an  arbi¬ 
trary  differential  pressure, 

ft  Vft?/sec 

Dynamic  pressure,  lc/ft^ 

Instantaneous  inflated  canopy 
surface  area,  ft" 

o 

Canopy  surface  area,  ft 


max 


-  Instantaneous  time,  sec 

-  Reference  time  when  the  parachute 
has  reached  the  design  drag  area 
for  the  first  time,  sec 

-  Canopy  inflation  time  ./hen  the 
Inflated  canopy  has  reached  its 
maximum  physical  size,  sec 

-  Air  velocity  through  cloth  in 
effective  porosity,  ft/sec 

-  Fictitious  theoretical  velocity 
used  in  effective  porosity,  ft/sec 

-  Instantaneous  system  velocity,  ft/sec 

-  System  velocity  at  the  time  t  =  t  , 

x't/sec  u 

-  System  velocity  at  the  end  of 
suspension  line  stretch,  ft/sec 

-  ‘tolune  of  air  which  must  be  collec¬ 
ted  during  the  inflation  process, 
ft^ 

-  Hardware  weight,  lb 

-  Instantaneous  shock  factor 


-  Shock  factor  at  the  time  t  «=  t 

-  Air  density,  slugs/ft’ 

-  patio  of  parachute  projected  mouth 
area  at  line  stretch  to  the  steady- 
state  projected  area 

-  Instantaneous  elongation 

-  Maximum  elongation 

-  Initial  elongation  at  the  beginning 
of  the  elastic  phase  of  Inflation 


S.F.  -  Para-hute  safety  factor  =  Fc/Fnax 
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Appendix  B 

A  GUIDE  FOR  THE  USE  OF  APPENDIX  A 


At  first  reading.  Appendix  A  may  appear  to  be  a  complicated  sys¬ 
tem  of  analysis  because  of  the  many  formulae  presented.  Actually, 
once  understood,  the  technique  is  straightforward  and  uncomplicated. 
The  author  has  attempted  to  simplify  the  algebra  wherever  possible. 
This  appendix  presents,  in  semi-outline  form,  a  guide  to  the 
sequence  of  calculations  because  the  analysis  does  require  use  of 
formulae  from  the  text,  not  necessarily  in  the  order  in  which  they 
were  presented.  Also,  the  user  can  be  referred  to  graphs  of  per¬ 
formance  to  illustrate  effects. 

In  order  to  compute  tQ,  other  parameters  must  be  obtained  from 
various  sources. 

I.  Determine  System  Parameters 

1.  CDS0,  drag  area,  ft2  obtained  from  design  requirement. 

2.  Vs,  fps,  velocity  of  system  at  suspension  line  stretch. 

3.  c,  slugs/ft3,  air  density  at  deployment  altitude. 

4.  W,  lb,  system  weight  (including  weight  of  the  parachute) 
from  design  requirements. 

5.  VQ,  ft3,  this  volume  of  air,  which  is  to  be  collected 
during  inflation,  is  calculated  from  the  steady-state  inflated 
shape  geometry  of  the  particular  parachute  type.  The  nomenclature 
is  described  in  Figure  24,p.A-14.  When  Dc  or  Dp  is  known,  a  can  be 
calculated  from  data  in  Table  I  and  Table  II,  p.  A-1S,  for  various 
parachute  types  and  number  of  gores.  Then  the  geometric  volume  V0 
can  be  calculated  by  Equation  (31),  p.  A-14,  with  appropriate 
values  of  b/a  and  b'/a  from  the  tables. 

6.  A^q,  ft2,  steady-state  canopy  mouth  area 
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where  N/a ,  b/a ,  and  b'/a  are  available  from  Tables  I  and  II  for  the 
particular  type  of  parachute  and  number  of  gores. 

7.  Asq,  ft2,  canopy  surface  area  = 

8.  C  ,  pressure  coefficient,  see  Figure  18,  p.  A-12.  A 
constant  Cp  =  1.7  for  all  altitudes  seems  to  yield  acceptable 
results . 


9.  Constants  k  and  n  are  derived  from  measurements  of  the  air 
flow  through  the  cloth.  Only  k  is  needed  for  Equation  (14),  but  n 
is  also  required  for  Equation  (13).  These  parameters  can  be  deter¬ 
mined  for  any  cloth  using  the  technique  described  beginning  on 
p.  A-12.  The  two-point  method  is  adequate  if  the  AP  across  the 
cloth  is  in  the  range  of  AP  for  actual  operation.  Check-points  of 
cloth  permeability  can  be  measured  and  compared  to  calculated 
values  to  verify  agreement.  If  the  data  are  to  be  extrapolated  to 
operational  AP ' s  greater  than  measured,  a  better  method  of  deter¬ 
mining  k  and  n  from  the  test  data  would  be  a  least  squares  fit 
through  many  data  points.  This  way  errors  due  to  reading  either 
of  the  two  points  are  minimized. 

11.  Step  1 

Calculate  the  reference  time  tQ  by  use  of  Equations  (13)  or 
(14),  p.  a-7.  If  the  deployment  altitude  is  50,000  feet  or  higher, 
Equation  (14)  is  preferred  due  to  its  simplicity.  For  altitudes 
from  sea  level  to  50,000  feet,  Equation  (13)  is  preferred.  Figure 

12,  p.  A-8,  shows  the  effect  of  altitude  on  tQ  and  can  be  taken  as 

a  guide  for  the  user  to  decide  whether  to  use  Equation  (13)  or  (14). 
One  should  keep  in  mind  that  the  opening  shock  force  can  be  a 
strong  function  of  inflation  time,  so  be  as  realistic  as  possible. 

If  Equation  (13)  is  elected,  the  method  in  use  at  the  NSWC/WO  is 
to  program  Equation  (13)  to  compute  the  parachute  volume,  Vq,  for 
an  assumed  value  of  tQ.  Equation  (14),  because  of  its  simplicity, 
can  be  used  for  a  first  estimate  of  tD  at  all  altitudes.  The  com- 
putc d  canopy  volume  is  then  compare!  to  the  canopy  volume  calculated 
from  the  geometry  of  the  parachute  as  per  Equation  (31),  p.  A-14. 

If  the  volume  computed  from  the  mass  flow  is  within  the  volume 
computed  from  the  geometry  within  plus  or  minus  a  specified  delta 
volume,  the  time  tG  is  printed  out.  If  not  within  the  specified 
limits,  t0  is  adjusted,  and  a  new  volume  calculated.  For  a  35-foot 
D0 ,  T-10  type  canopy,  I  use  plus  or  minus  10  cubic  feet  in  the 
volume  comparison.  The  limit  would  be  reduced  for  a  parachute  of 


sma  Her 

Do 

• 

If 

Vo 

calculated 

=  v0 

geometry  ± 

10, 

then 

print  answer. 

If 

v0 

calculated 

*  VD 

geometry  + 

10, 

then 

correct  tQ  as  follows 

Y0  9#om*try 

*o  y  calculated  (B-2) 

O 
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The  new  value  of  tQ  is  substituted  in  the  'do  loop”  and  the 
volume  recomputed.  This  calculation  continues  until  the  required 
volume  is  within  the  specified  limits. 

HI.  Calculate  tQ  corrected  for  initial  area.  The  tQ  of  Section 
II  assumes  that  the  parachute  inflated  from  a  zero  initial  area. 

If  this  is  a  reasonable  assumption  for  the  particular  system  under 
study,  then  the  mass  ratio  can  be  determined  from  Equation  (6), 
p.  A- 4 .  For  n  =  0  if  the  value  of  M  i.  0.19,  then  a  finite  state 
of  deployment  exists,  and  the  time  ratio  of  occurrence  and  the 
maximum  shock  factor  can  be  determined  from  Equations  (9)  and  (10), 
respectively,  on  p.  A-5.  If  n  4  0,  then  the  limiting  mass  ratio 
for  finite  operations  will  rise  slightly  as  described  in  Appendix 
C.  Figures  C-l  and  C-2  illustrate  the  effects  of  initial  area  on 
limiting  mass  ratios  and  shock  factors  respectively.  If  the  mass 
ratio  is  greater  th  a  n  the  limiting  mass  ratio  (ML) ,  then  the 
maximum  shock  force  occurs  at  a  time  greater  than  tD  and  the 
elasticity  of  the  materials  must  be  considered  (see  Section  VI). 

If  n  /  0,  then  the  reference  time,  t  ,  will  be  reduced,  and 
the  mass  ratio  will  rise  due  to  partial  inflation  at  the  line 
stretch.  Figures  9  and  10,  p.  A-6,  illustrate  the  effects  of 
initial  area  on  the  velocities  and  shock  factor  during  the 
"unfolding”  inflation.  Equation  (15),  p.  A-9,  can  be  used  to 
correct  tQ  calculated  for  the  cases  where  n  =  Ai/Ac.  If  the  initial 
value  of  drag  area  is  known.  Equation  (16)  ,  p.  A-9,  can  be  used  to 
correct  t  and  rechecked  for  limiting  mass  ratios  versus  n  in 
Appendix  C. 

IV.  Opening  shock  calculations  in  the  elastic  phase  of  inflation. 

It  has  been  considered  that  from  time  t  =  0  to  t  =  tQ  the  para¬ 
chute  has  been  inelastic.  At  the  time  t  =  t0  the  applied  aero¬ 
dynamic  load  causes  the  materials  to  stretch  and  the  parachute 
canopy  increases  in  size.  The  increased  size  results  in  an 
increase  in  load,  which  causes  further  growth,  etc.  This  sequence 
of  events  continues  until  the  applied  forces  have  been  balanced  by 
the  strength  of  materials.  The  designer  must  insure  that  the 
constructed  strength  of  the  materials  is  sufficient  to  resist  the 
applied  loads  for  the  material  elongation  expected.  Use  of 
materials  of  low  elongation  should  result  in  lower  opening  shock 
forces  as  CDSmax  is  reduced. 

When  the  mass  ratio  of  the  system  is  greater  than  the  limiting 
mass  ratio,  the  elasticity  of  the  materials  and  material  strength 
determine  the  maximum  opening  shock  force.  The  maximum  elongation 
£max  an<^  t*ie  ultimate  strength  of  the  materials  are  known  from 
tests  or  specifications.  The  technique  begins  on  p.  A-9. 

At  the  time  t  =  tc,  calculate  the  following  quantities  for  the 
particular  values  of  M  and  n. 

a.  V0/Vs  from  Equation  (18),  p.  A-9. 
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b.  XQ  from  Equation  (22),  p.  A-10. 

c.  eQ  from  Equation  (23),  p.  A-10. 

d.  Determine  CDSmax/CDS0  from  Figure  15,  p.  A-10. 

e.  Calculate  the  inflation  time  ratio  tf/t~  from  Equation 

(24) ,  p.  A-10.  ° 

f.  Calculate  the  maximum  shock  factor  from  Equation  (19), 
p.  A- 9. 

g.  Calculate  the  opening  shock  force  Fm«.  “  where 


F»  ■  i»'Scos. 


h.  Calculate  filling  time,  tf(sec) 


V.  In  order  to  simplify  the  required  effort,  the  work  sheets  of  Table 
B-l  are  included  on  pages  B-5  through  B-9  to  aid  the  engineer  in 
systematizing  the  analysis.  The  work  sheets  should  be  reproduced  to 
provide  additional  copies. 
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Table  B-l  .  Op>ening  Shock  Force 

CALCULATION  WORK  SHEETS 


1.  Parachute  type  - 

2.  System  parameters 

a.  System  weight,  W  (lb) 

2 

b.  Gravity,  g  (ft/sec  ) 

c.  Deployment  altitude  (ft) 

d.  Deployment  air  density,  C  (slugs/ft3) 

e.  Velocity  at  line  stretch,  Vg  (fps) 

f.  Steady  state  canopy  data 


(1) 

Diameter,  D0 

(ft) 

(2) 

Inflated  diameter,  2a  (ft); 

O) 

Surrace  area 

rr  _  2 

,  S0  (ft-);  -  "o 

(4) 

Drag  area,  C 

Ds0  (ft2);  CDx 

(5) 

Mouth  area, 

am*0  (ft2) 

Amo  '  »' J’  * 

(6) 

Volume,  V,,* 

(ft3) 

*'+•**  [*•*] 


p 

V. 


2i 


A,o‘So 

C„S 
D  o 


amo 


lb. 

ft/»c2 

ft. 

Ip*. 

ft. 

ft. 

ft.2 

ft.2 


ft. 2 


ft.J 


g.  Cloth  data 

(1)  k  )  Calculate  using  technique  beginning  on 

(2)  n  )  p.  A-12.  Note:  Permeab il ity  is  usually 

measured  as  f t 3/ f l ^ /min . 

For  these  calculations 
permeability  must  be 
expressed  as  ft-Vft^/sec 


k 

n 


*  Data  for  these  calculations  are  listed  in  Tables  1 
and  2  ,  p.  A- 1 5 . 


B  6 
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Table  B-l .  Opening  Shock  Force 

(cont 'd) 

(3)  £max'  determine  maximum  elongations  from 

pull  test  data  of  joints,  seams,  lines,  etc.  Use 
minimum  determined  from  tests. 

(4)  Cp;  pressure  coefficient 

h.  Steady  state  drag,  Fs  (lb),  Fg  *  2pVSC0So 

i.  Parachute  constructed  strength,  Fc  (lb);  deter¬ 
mined  from  data  on  efficiency  of  seams,  joints,  lines. 
Constructed  strength  is  the  minimum  load  required  to 
fail  a  member  times  the  number  of  members. 

3.  Force  calculations 

a.  Calculate  tQ  for  n  *  0;  eq.  14,  p.  A-7. 


r  P*Y°  r 

c°*«  In 

14W 

t  m 

2W  L  Amo  -  Ago  k  /cPp\  i  J 

°  P«VSC0So 

L< 

Vt/  -ij 

Check  Figure  13, 
using  eq.  13,  p.  A-7. 

p.  A-8,  for 

advisability  of 

b.  If  n  -  0,  proceed  with  steps 
If  n  f  0,  go  to  step  f. 

c  through  e. 

c.  Mass  ratio,  M;  eq.  6,  p.  A-4 

M 


2W 


POW  os. 


d.  If  M  <  4/21  for  n  =  0,  then  finite  mass 
deployment  is  indicated. 

(1)  Time  ratio  at  x^  eq>  9,  p-  A- 3 


2 

-L  .  (mf 

‘.♦x,  \  4  I 


(2)  Max  shock  factor,  eq. ,  10,  p.  A- 5 


X, 


max 


16  (  _21m\* 

49  V  4  ) 


B-6 


'-P 

F. 


M 


t 

®Xj 


mtx 


max 


lb. 


lb. 


NSWC  TR  87  96 


(3) 


Table  B-l 
Max  .shock  force, 


Opening  Shock 
Fmax  <lb> 


Force 

(Cont'd) 


lb. 


e.  If  M  >  4/21;  then  Intermediate  mass  or  infinite 
mass  deployment  is  indicated  and  the  elasticity  of 
materials  is  involved.  Calculate  the  trajectory  con¬ 
ditions  at  time  t  =  tQ. 


(1)  Velocity  ratio  @  t  =  tQ  for  n  =  0 


V 

o 


1 


+ 


1 

7M 


V 


(2)  Shock  factor  XQ  Q  t  =  tQ  for  n  =  0 


*o  - 


['•*] 


2  I  V„ 


(3)  Initial  elongation,  CQ;  eq.  23,  p.  A-10 


xVs 


p.  A-LO 


cns. 


(4)  Determine  ,va—  from  Figure  15,  p.  A-  10 


"Do 


(5)  Calculate  inflation  time  ratio,  tf  ;  eq .  24, 


CDStr.»n 


co% 


19,  p.  A- 9 


(6)  Calculate  maximum  shock  factor,  raax;  eq . 


/  vy 

Vo/ 


""  K-MW-ir 

(7)  Calculate  maximum  shock  force,  I  .  '"A.-) 


x 


o 


c 


o 


cDsmax 

CDSo 


*o 


max 


lb. 
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Table  B-l .  Opening  Shock  Force 

(cont ' d ) 


(8)  Inflation  time,  sec  =  “  *o 


(S) 


f.  If  n  t  0,  correct  tn  for  initial  area  effects; 
eq.  16,  p.  A- 9 


g.  Mass  Ratio,  M,  eq.  6,  p.  A-4 

2W 


■-&n 


t0  calculated 


M 


P9Vs«„CdSo 


h.  Calculate  limiting  mass  ratio, 

1 


M 


L  3(1  -  rj) 


9  3  1 

U"  14  72  7 


If  M  <  ,  finite  mass  deployment  is  indicated  and  x^ 

can  be  determined  by  eq.  8,  p.  A-3  by  assuming  values 
of  t/tQ  and  plotting  the  data  using  the  methods  of 
Appendix  C. 


i.  If  M  >  Ml,  then  intermediate  mass  or 
infinite  mass  deployment  is  indicated  and  the 
elasticity  of  materials  is  involved.  Calculate 
the  trajectory  conditions  at  time  t  =  t0. 


p.  A- 9 


(1)  Velocity  ratio  (?  t  =  tQ  for  n  5*  0;  eq .  18, 


V„ 


22,  p .  A- 1 0 


*0  - 


r,.i 

r<.L-II>!.  +  ullznl  +  r,2 

L  » 

L  7  2  ^  J 

factor  XQ  f?  t  -  tQ  for 

1 

1*1 

IL-sii  *  — il 

M 

7  2 

(3)  Initial  elongation,  f  ;  eq.  23,  p.  A-10 


Vs 


V  “ 


Cr,S, 


(4)  Determine  ■  —  ^,aj<  from  Figure  15,  p.  A-10 


CDSo 


B  8 


X 

tf 

*0 

M 

M, 

L 

Vo 

v, 

*0 

£ 

0 

C0SfP8X 

cDs0 

S«c. 


S«c. 
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Table  B-l.  Opening  Shock  Force  (ConbasA 


(5)  Calculate  inflation  tine  ratio,  —  ;  eq 
p.  A- 10  ,  r° 


>.  \  v. 


(6)  Calculate  maximum  shock  factor,  max'» 
eq.  19,  p.  a-9  / t  \ 6 


(7)  Calculate  maximum  shock  force,  F  y(lb) 


F  »  x  Fc 
max  vi_  S 
max 


(8)  Calculate  inflation  time,  tf(sec) 


V  *  «o 
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The  calculated  canopy  volume,  Us  calc»  is  determined  from  equation 

(11). 


h  •  X  [vA”°(t)6  -  A,<0 k  v2"] 


A  program  for  calculating  tn  for  solid  cloth  parachutes  and 
the  opening  shock  force  profile  during  the  inflation  of  several 
parachute  types  is  provided  in  Table  1.  Equations  (4),  (11),  and 
(12)  are  programmed  together  with  the  vertical  deployment  opening 
shock  equations  (2)  through  (8)  in  FORTRAN  IV  language.  The 
included  examples  were  calculated  via  the  program  using  a  VAX  780 
computer.  The  program  operates  in  two  modes.  Mode  1,  for  solid 
cloth  parachutes,  calculates  the  vertical  deployment  reference 
time  t0  for  the  parachute  system  parameters  and  operational 
deployment  data,  and  then  calculates  the  opening  shock  profile 
during  inflation.  A  typical  data  print  out  is  shown  in  Table  2. 

It  is  necessary  to  estimate  an  initial  value  of  "t0".  The  program 
calculates  the  canopy  volume  for  the  estimated  time  and  compares 
the  Up  calc  to  the  volume  derived  from  the  canopy  geometry.  If  the 
calculated  volume  is  not  within  specified  limits",  the  program 
adjusts  “t0"  by  equation  (12)  and  reiterates  the  program  until  the 
calculated  volume  is  within  the  specified  limits. 


Mode  2  of  the  program  calculates  opening  shock  profiles  for  input 
values  of  tQ.  Mode  2  analysis  of  other  types  of  parachutes  is 
possible  by  the  selection  of  the  proper  values  of  MjH  (1/2,  1,  2, 
3,  4,  5,  or  6)  and  **tm.  The  opening  shock  force  variation  for 
examples  (1)  and  (2)  are  plotted  in  Figures  3  and  4.  The  nominal 
t0  for  n«0.632  was  calculated  by  the  program  in  mode  1  and  the 
force-time  survey  was  calculated  in  mode  2. 

Figure  3  illustrates  that  j  jchutes  deployed  in  a  vertical 
toward-the-earth  trajectory  inflate  faster  than  the  same  system 
deployed  horizontally  at  the  same  altitude  and ^velocity. 

Inflation  reference  times  for  parachute  types  other  than 
solid  cloth  canopies  can  be  developed  from  the  mass  flow  equation. 
Thi6  requires  that  the  flow  through  the  canopy  be  expressed  in  a 
form  similar  to  the  solid  cloth  canopy  cloth  permeability,  P, 
where  the  rate  of  flow  per  unit  area  is  a  function  of  the  pressure 
differential  across  the  cloth  or  grid. 


13 


Change  1 


NSWC  TR  87-96 


TABLE  1.  INSTANTANEOUS  DRAG  AREA,  VELOCITY.  OPENING  SHOCK  FORCE.  AND 
DISTANCE  OF  FALL  OF  A  PARACHUTE  DEPLOYED  IN  VERTICAL  FALL. 


THIS  PROGRAM  CALCULATES  THE  INSTANTANEOUS  DRAG  AREA,  VELOCITY, 
OPENINO  SHOCK  FORCE,  AND  DISTANCE  OF  FALL  OF  A  PARACHUTE  DEPLOYED  IN 
VERTICAL  FALL 


THE  PROGRAM  OPERATES  IN  TWO  MOOES: 

MODE  1  -  CALCULATES  THE  INFLATION  TIME  AND  PERFORMANCE  PROFILES  FOR  SOLID  CLOTH 

PARACHUTES  (TO  INPUT  AS  INITIAL  ESTIMATE)  (IOPT  =  1) 

MODE  2  -  CALCULATES  THE  PERFORMANCE  PROFILES  FOR  VARIOUS  TYPES  OF  PARACHUTES  (J) 

INFLATION  TIME  INPUT  IS  REQUIRED  (IOPT  =  2) 


INPUT  IOPT  -  1  <FQR  MODE  1) 

-  2  (FOR  MODE  2) 

INPUT  NEEDED  FOR  BOTH  MODES 

RHO  -  AIR  DENSITY  AT  GIVEN  ALTITUDE  (SLUGS/FT3) 

VS  -  VELOCITY  AT  SUSPENSION  LINE  STRETCH  (FT/SEC) 

CDSO  -  DESIGN  DRAG  AREA  (FT2) 

TO  -  I OPT* l  INITIAL  GUESS  FOR  INFLATION  REF  TIME  (SEC) 

I 0PT«2  ACTUAL  INFLATION  REFERENCE  TIME  (SEC) 

U  -  WEIGHT  (LUS) 

J  -  *6  FOR  FLAT  CIRCULAR  PARACHUTE 

-1  FOR  RIBBON  TYPE  OF  PARACHUTE 


/); 

1  I 


{ 

1 

\ 

« 

A 


INPUT  NEEDED  FOR  IOPT  =  1  ONLY 


AMO  -  STEADY-STATE  MOUTH  AREA  ( FT 2 ) 

ASO  -  CANOPY  DESIGN  SURFACE  AREA  ( FT2  > 
K  -  CLOTH  PERMEABILITY  CONSTANT 

CP  -  PRESSURE  COEFFICIENT 

N  -  CLOTH  PERMEABILITY  EXPONENT 

VO  -  GEOMETRIC  VOLUME  (FT3 ) 


REAL*4  N 
TODEN* lOOOOO 
5  PRINT  *,  ' INPUT  IOPT  ' 

READ <3,  *,  END*- 1 00!  IOPT 

PRINT  *.  'INPUT  RHO, VS, CDSO, TO, W, J' 

READ( 3,  * ) RHO,  VS,  CDSO.  TO,  W,  J 

IF ( IOPT  EQ  2)  CO  TO  3 

PRINT  *,  ' INPUT  AMO. ASO.  XK. CP,  N,  VO ' 

READ (3,  * ) AMO. ASO, XK,  CP,  N, VO 
3  DT-TO/TODEN 

TAU-0  j 

0‘32  2 

S=0  | 

X*TAU 

CDS»TAU*CDSO 

FS-  3«RH0*VS**2*CD50  j 

F«TAU*FS 

VOL-O 

IPA5S*0 

l 

/ 


1 

!  J 


1A 


